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Measurement of the effective 



This paper reports the first measurement of the effective B® — > J/ip K® lifetime and 
an updated measurement of its time-integrated branching fraction. Both measure- 
ments are performed with a data sample, corresponding to an integrated luminosity 
of 1.0 fb _1 of pp collisions, recorded by the LHCb experiment in 2011 at a centre- 
of-mass energy of 7 TeV. The results are: Tj% K o = 1.75 ± 0.12 (stat) ± 0.07 (syst) ps 

and B(B° S -)• JfyK®) = (1.97 ± 0.23) x 10~ 5 . For the latter measurement, the un- 
certainty includes both statistical and systematic sources. 
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1 Introduction 



In the Standard Model (SM), CP violation arises through a single phase in the CKM 
quark mixing matrix [Tj . In decays of neutral B mesons (B stands for a B° or B® meson) 
to a final state accessible to both B and B, the interference between the amplitude for 
the direct decay and the amplitude for decay via oscillation leads to time-dependent CP 
violation. A measurement of the time-dependent CP asymmetry in the B° — > J/ipK® 
mode allows for a determination of the B°-B° mixing phase <pd- In the SM it is equal 
to 2/3 |2), where is one of the angles of the unitarity triangle in the quark mixing 
matrix. This phase has already been well measured by the B factories |3}|4], but further 
improvements are still necessary to conclusively resolve possible small tensions with the 
other measurements constraining the unitarity triangle pi. The latest average composed 
by the Heavy Flavour Averaging Group (HFAG) is sin fa = 0.682 ± 0.019 (6|. To achieve 
precision below the percent level, knowledge of the doubly Cabibbo-suppressed higher 
order perturbative corrections, originating from penguin topologies, becomes mandatory. 
These contributions are difficult to calculate reliably and therefore need to be determined 
directly from experimentally accessible observables. 

From a theoretical perspective, the B° s — > J/ip K® mode is the most promising candidate 
for this task. It is related to the B° — > J/ipK® mode through the interchange of all 
d and s quarks ([/-spin symmetry, a subgroup of SU(3)) [7], leading to a one-to-one 
correspondence between all decay topologies of these two modes, as illustrated in Fig. [T] 
Moreover, the B° s ->• J/i/jK° penguin topologies are not CKM suppressed relative to the 
tree diagram, as is the case for their B° counterparts. A further discussion regarding the 
theory of this decay and its potential use in LHCb is given in Ref. (8l. 

To determine the parameters related to the penguin contributions in these decays, a 
time-dependent CP violation study of the B° s — > J/ipK® mode is required. The determi- 
nation of its branching fraction, previously measured by CDF [9] and LHCb [To] , was an 
important first step, allowing a test of the [/-spin symmetry assumption that lies at the 
basis of the proposed approach. The second step towards the time-dependent CP vio- 
lation study is the measurement of the effective B° s — > J/ip lifetime, formally defined 
as 11 

eff _ fZ°t(T(B a (t)-+Jftl,K°))tit 
Tj#k%- f o °°(T(B s (t)^J/i/;K$))dt ' [L) 

where 

(T(B s (t) -+ JftK*)) = T(B s (t)^J/ijK s ) + T(B° s (t)^J/^K°) (2) 

= R K e- rHt + R L e~ Tht (3) 

is the untagged decay time distribution, under the assumption that CP viola- 
tion in B° s -B° s mixing can be neglected 6|. Due to the non-zero decay width 
difference Ar s = T H - T L = 0.106 ± 0.013 ps^ (l2) between the heavy and light B° s 
mass eigenstates, the effective lifetime does not coincide with the B® lifetime 



tbo = l/r s = 1.513 ± 0.011 ps 12 , where T s = (Tn + Tl)/2 is the average B® decay 



1 




width. Instead, it depends on the decay mode specific relative contributions i?n and R^. 
These two parameters also define the CP observable 



A 



AF 3 



Ru — Rh 

Rh + Rh 



which allows the effective lifetime to be expressed as 11 

1 + 2 A AFs y s + vl 



eff 



TBO 



yi 



1 + -AAr 3 ys 



(4) 



(5) 



where y s = Ar s /2r s is the normalised decay width difference. For the B® — > J/ipK® 
mode, the value of ^4.Ar s depends on the penguin contributions, and in particular on their 

. Using the latest estimates on the size of the B° s — > J/ipK® 
gives -4 A r s = 0.944 ± 0.066 and the SM prediction 



relative weak phase cj) s ff 
penguin contributions 13 



eff 



SM 



1.639 ±0.022 ps 



(6) 



Effective lifetime measurements have been performed for the B° s — > K + K 14 and 
B° s J/^/ (980) 1 15] decay modes. 

This paper presents the first measurement of the effective B° s — > J/ipK® lifetime, as 
well as an update of the time-integrated branching fraction measurement in Ref. [lO] , 
performed with a data sample, corresponding to an integrated luminosity of 1.0 fb~ of 
pp collisions, recorded at a centre-of-mass energy of 7 TeV by the LHCb experiment in 
2011. 

The LHCb detector [16] is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < r\ < 5, designed for the study of particles containing b or c quarks. 
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The detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4Tm, and three stations 
of silicon-strip detectors and straw drift tubes placed downstream. The combined track- 
ing system has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% 
at 100 GeV/ c, and impact parameter resolution of 20 \xm for tracks with high transverse 
momentum (pp) with respect to the beam direction. Charged hadrons are identified using 
two ring-imaging Cherenkov detectors [17]. Photon, electron and hadron candidates are 
identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 
an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 
system composed of alternating layers of iron and multiwire proportional chambers. 

Events are selected by a trigger system 18 consisting of a hardware trigger, which 
requires muon or hadron candidates with high p^, followed by a two-stage software trigger. 
In the first stage a partial event reconstruction is performed. For this analysis, events 
are required to have either two oppositely charged muons with combined mass above 
2.7 GeV/c 2 , or at least one muon or one high-pr track (p<i > 1.8 GeV/c) with a large 
impact parameter with respect to all pp interaction vertices (PVs). In the second stage a 
full event reconstruction is performed and only events containing J/ip — > /i + fjT candidates 
are retained. 

The signal simulation samples used for this analysis are generated using 



Pythia 6.4 19 with a specific LHCb configuration 20 . Decays of hadronic particles are 



described by EvtGen [21 1 in which final state radiation is generated using Photos 22 



The interaction of the generated particles with the detector and its response are imple- 



mented using the Geant4 toolkit 23 as described in Ref. 24 



2 Data samples and initial selection 

Candidate B — > J/ipK® decays are reconstructed in the J/ip — > and K® — > 7r + 7r~ 

final state. Candidate J/ip — > /x + \T decays are required to form a good quality vertex and 
have a mass in the range [3030,3150] MeV/c 2 . This interval corresponds to about eight 
times the mass resolution at the J/ip mass and covers part of the J/ip radiative tail. 

The selected J/ip candidate is required to satisfy the trigger decision at both software 
trigger stages. The K® selection requires two oppositely charged particles reconstructed 
in the tracking stations placed on either side of the magnet, both with hits in the vertex 
detector ('long candidate) or without ('downstream K®' candidate). The long (down- 
stream) — > 7r + 7r~ candidates are required to form a good quality vertex and have a 



■■s 

mass within 35 (64) MeV/c 2 of the known K® mass 25 . Moreover, to remove contami- 
nation from A pn~ decays, the reconstructed pir~ mass of the long (downstream) K® 
candidates is required to be more than 6 (10) MeV/c 2 away from the known A mass [25] . 
Furthermore, the K® candidates are required to have a flight distance that is at least five 
times larger than its uncertainty. 

Candidate B mesons are selected from combinations of J/ip and K® candidates with 
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mass ifij^j-o in the range [5180, 5520] MeV/c 2 . The reconstructed mass and decay time are 
obtained from a kinematic fit [26) that constrains the masses of the and 7r + 7r~ pairs 

to the known J/0 and K® masses [25], respectively, and constrains the B candidate to 
originate from the PV. In case the event has multiple PVs, all combinations are considered. 
The x 2 °f the fit, which has eight degrees of freedom, is required to be less than 72 and the 
estimated uncertainty on the B mass must not exceed 30 MeV/c 2 . Candidates are required 
to have a decay time larger than 0.2 ps. To remove misreconstructed B° — > J/ipK*° 
background that survives the requirement on the K® flight distance, the mass of the long 
B° — > J/ipK® candidates computed under the J/ipK ± n T mass hypotheses must not be 



within 20 MeV/c 2 of the known B° mass 25 



3 Multivariate selection 

The loose selection described above does not suppress the combinatorial background suffi- 
ciently to isolate the small B° s — > J/ipK® signal. The initial selection is therefore followed 
by a multivariate analysis, based on a neural network (NN) [27]. The NN classifier's 
output is used as the final selection variable. 

The NN is trained entirely on data, using the B° — > J/ipK® signal as a proxy 
for the B® — > J/ipK® decay. The training sample is taken from the mass windows 
[5180, 5340] MeV/c 2 and [5390, 5520] MeV/c 2 , thus avoiding the B° s signal region. A nor- 
malisation sample consisting of one quarter of the candidates, selected at random, is left 
out of the NN training to allow an unbiased measurement of the B° yield. The signal 
and background weights are determined using the sPlot technique [28] and obtained by 
performing an unbinned maximum likelihood fit to the mass distribution of the candi- 
dates surviving the loose selection criteria. The fitted probability density function (PDF) 
is defined as the sum of a B° signal component and a combinatorial background. The 
parametrisation of the individual components is described in more detail in the next 
section. 

Due to the differences in the distributions of the input variables of the NN, as well 
as the different initial signal to background ratio, the multivariate selection is performed 
separately for the B candidate samples containing long and downstream K® candidates. 
In the remainder of this paper, these two datasets will be referred to as the long and down- 
stream K® sample, respectively. The NN classifiers use information about the candidate 
kinematics, vertex and track quality, impact parameter, particle identification informa- 
tion from the RICH and muon detectors, as well as global event properties like track and 
interaction vertex multiplicities. The variables that are used in the NN are chosen to 
avoid correlations with the reconstructed B mass. 

Final selection requirements on the NN classifier outputs are chosen to optimise the 
expected sensitivity to the B° s signal observation. The expected signal and background 



yields entering the calculation of the figure of merit 29 are obtained from the normalisa- 
tion sample by scaling the number of fitted B° candidates, and by counting the number 
of events in the mass ranges [5180,5240] MeV/c 2 and [5400,5520] MeV/c 2 , respectively. 
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Long K® Downstream K® 




Figure 2: Fitted B — > J/ipK® candidate mass distributions and their associated residual 
uncertainties (pulls) for the (left) long and (right) downstream K® samples, after applying 
the final requirement on the NN classifier outputs. 



After applying the final requirement on the NN classifier output associated with the long 
(downstream) K® sample, the multivariate selection rejects, relative to the initial selec- 
tion, 98.7% (99.6%) of the background while keeping 71.5% (50.2%) of the B° signal. Due 
to the worse initial signal to background ratio, the final requirement on the NN classifier 
output is much tighter in the downstream K® sample than in the long K® sample. 

After applying the full selection, the B candidate can still be associated with more 
than one PV in about 1% of the events. Likewise, about 0.1% of the selected events have 
several candidates sharing one or more tracks. In these cases, respectively one of the 
surviving PVs and one of the candidates is used at random. 



4 Event yields 

For the candidates passing the NN requirements, the ratio of B® and B° yields is deter- 
mined from an unbinned maximum likelihood fit to the mass distribution of the recon- 
structed B candidates. The fitted PDF is defined as the sum of a B° signal component, a 
B Q S signal component and a combinatorial background. The B® component is constrained 



to have the same shape as the B° PDF, shifted by the known B®-B° mass difference 30 
The mass lineshapes of the B — > J/ipK® modes in both data and simulation exhibit 
non-Gaussian tails on both sides of their signal peaks due to final state radiation, the 
detector resolution and its dependence on the decay angles. Each individual signal shape 



is parametrised by a double-sided Crystal Ball (CB) function 31 . The parameters de- 
scribing the CB tails are taken from simulation; all other parameters are allowed to vary 
in the fit. The background contribution is described by an exponential function. 
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Table I: Signal yields from the unbinned maximum likelihood fits to the B -> J/i[>K° 
candidate mass distributions. The uncertainties are statistical only. The yield ratio is 
calculated from the quantities highlighted in boldface, where the fitted B° yield is first 
multiplied by a factor of four. 



Sample 


Yield 


Long Kg 


Downstream K® 


Normalisation 


B° J/il>K® 

B° s -> J/i/jK° s 
Background 


2205 ± 47 

21 ± 5 
56 ± 11 


3651 ± 61 

49 ± 8 
110 ± 16 


Full 


B° JfaKl 
B° s -> J/ij>K° 
Background 


9031 ± 96 
115 ± 12 

287 ± 23 


14,391 ± 122 
158 ± 15 

490 ± 32 


Yield ratio R = 


ArPull I a /\T Norm 
B®—>-J/4:Kgi B —>-J/'4>Kg 


0.0131 ± 0.0014 


0.0108 ± 0.0010 



Average yield ratio R 0.0116 ± 0.0008 



The results of the fits are shown in Fig. [2j and the fitted yields are listed in Table [T] 
The B° yield is determined in the normalisation sample and scaled to the full sample, 
whereas the B° s yield is obtained directly from the full sample. The scaled B° yield, 
obtained from the unbiased sample, differs from the corresponding fit result in the full 
sample by —211 ± 211 events for the long sample and by 213 ± 273 events for the 
downstream K® sample. Both results are in good agreement, showing that the NN is not 
overtrained. The yield ratios obtained from the long and downstream K® samples are 
compatible with each other and are combined using a weighted average. 



5 Decay time distribution 



Following the procedure explained in Ref. 32 , the effective B® — > J/ip K® lifetime is 



determined by fitting a single exponential function g{t) oc exp(— t/r s i ng i e ) to the decay time 
distribution of the B® — > J/ipK® signal candidates. In this analysis, the exponential shape 
parameter r sing i e is determined from a two-dimensional unbinned maximum likelihood fit 
to the mass and decay time distribution of the reconstructed B candidates. The fitted 
PDF is again defined as the sum of a B° signal component, a B° s signal component and 
a combinatorial background. The freely varying parameters in the fit are the signal and 
background yields, and the parameters describing the acceptance, mass and background 
decay time distributions. 

The decay time distribution of each of the two signal components needs to be corrected 
with a decay time resolution and acceptance model to account for detector effects. The 
shape of the acceptance function affecting the B° s — > J/ipK® mode is, like the lineshape 
of its mass distribution, assumed to be identical to that of the B° — > J/ipK® component. 
The acceptance function is obtained directly from the data using the B° — > J/tpK® mode. 
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B decay time [ps] B decay time [ps] 



Figure 3: Fitted B — > J/ipK® candidate decay time distributions and their associated 
residual uncertainties (pulls) for the (left) long and (right) downstream samples, after 
applying the final requirement on the NN classifier outputs. 



Contrary to the B® system, the B° system has a negligible decay width difference Ar^ [25] . 
The decay time distribution of the B° — > J/ip K® channel is therefore fully described by a 
single exponential function with known lifetime r^o = 1.519 ps [H]. Hence, fixing the B° 
lifetime to its known value allows the acceptance parameters to be determined from the 
fit. 

From simulation studies it is found that the decay time acceptance of both signal 
components is well modelled by the function 

Ace (t) = • (7) 



The parameter (3 describes the fall in the acceptance at large decay times 12 . The param- 
eters k and A model the turn-on curve, caused by the use of decay time biasing triggers, 
the initial selection requirements and, most importantly, the NN classifier outputs. 

The decay time resolution for the signal and background components is determined 
from candidates that have an unphysical, negative decay time. Due to the requirement of 
0.2 ps on the decay time of the B candidates applied in the initial selection, such events 
are not present in the analysed data sample. Instead, a second sample, that is prescaled 
and does not have the decay time requirement, is used. This sample consists primarily of 
J/ip mesons produced at the PV which are combined with a random K® candidate. The 
decay time distribution for these events is a good measure of the decay time resolution 
and is modelled by the sum of three Gaussian functions sharing a common mean. Two 
of the Gaussian functions parametrise the inner core of the resolution function, while the 
third describes the small fraction of outliers. 
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Figure 4: Fitted B — > J/ipK® candidate decay time distributions and their associated 
residual uncertainties (pulls) for the (left) long and (right) downstream candidates in 



the B® signal region mj^ K o 
on the NN classifier outputs. 



G [5340,5390] MeV/c 2 , after applying the final requirement 



The background decay time distributions are studied directly using the data. Their 
shape is obtained from background candidates that are isolated using the background 
weights determined by the sPlot technique, and cross-checked using the high mass side- 
band. The exact values of the shape parameters are determined in the nominal fit. Be- 
cause of the differences induced by the multivariate selection, the background decay time 
distribution of the long and downstream Kg samples cannot be parametrised using the 
same background model. For the long K® sample, the background is modelled by two 
exponential functions, describing a short-lived and a long-lived component, respectively. 
In the downstream K® sample such a short-lived component is not present due to the 
tighter requirement on the NN classifier output. Its decay time distribution is better de- 
scribed by a single exponential shape corrected by the acceptance function in Eq. (|7p with 
independent parameters («', a', /?'). The parameter f3' is set to zero because we also fit the 
lifetime of the single exponential function itself, and the combination of both parameters 
would result in ambiguous solutions. 

The decay time distributions resulting from the two-dimensional fits are shown in 
Figs. |3] and |4] for candidates in the full mass range mj/^ K o G [5180,5520] MeV/c 2 and 
in the B® signal region mj/^ K o G [5340,5390] MeV/c 2 , respectively. The fitted values 
are r sing i e = 1.54 ± 0.17 ps and r single = 1.96 ± 0.17 ps for the long and downstream 
Kg sample, respectively. The 1.7a difference between both results is understood as a 
statistical fluctuation. The two main fit results are therefore combined using a weighted 
average, leading to 

Single = 1.75 ± 0.12 ps , 
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where the uncertainty is statistical only. The event yields obtained from the two- 
dimensional fits are compatible with the results quoted in Table [TJ 

6 Corrections and systematic uncertainties 

A number of systematic uncertainties affecting the relative branching fraction 
B(B° S -»■ J/i/)K°)/B(B -»■ J/^Kl) and the effective lifetime are considered. The 
sources affecting the ratio of branching fractions are discussed first, followed by those 
contributing to the effective lifetime measurement. 

The largest systematic uncertainty on the yield ratio comes from the mass shape model, 
and in particular from the uncertainty on the fraction of the B° — > J/ipK® component's 
high mass tail extending below the B° s signal. The magnitude of this effect is studied by 
allowing both tails of the CB shapes to vary in the fit. The largest observed deviation in 
the yield ratios is 3.4%, which is taken as a systematic uncertainty. The mass resolution, 
and hence the widths of the CB shapes, is assumed to be identical for the B° and B® 
signal modes, but could in principle depend on the mass of the reconstructed B candidate. 
This effect is studied by multiplying the CB widths of the B® signal PDF by different 
scale factors, obtained by comparing B° and B® signal shapes in simulation. The largest 
observed difference in the yield ratios is 1.4%, which is taken as a systematic uncertainty. 
Varying the B° s -B° mass difference within its uncertainty has negligible effect on the yield 
ratios. 

The selection procedure is designed to be independent of the reconstructed B mass. 
Simulated data is used to check this assumption, and to evaluate the difference in selection 
efficiency arising from the different shapes of the B° — > J/tpK® and B° s — > J/tpK® decay 
time distributions. The ratio of total selection efficiencies is equal to 0.968 ± 0.007, and 
is used to correct the yield ratio. 

The stability of the multivariate selection is verified by comparing different training 
schemes and optimisation procedures, as well as by calculating the yield ratios for different 
subsets of the long and downstream K® sample. All of these tests give results that are 
compatible with the measured ratio. 

The corrections and systematic uncertainties affecting the branching fraction ratio 
are listed in Table [2] The total systematic uncertainty is obtained by adding all the 
uncertainties in quadrature. 

The main systematic uncertainties affecting the effective B° s — > J/ipK® lifetime arise 
from modelling the different components of the decay time distribution. Their ampli- 
tudes are evaluated by comparing the results from the nominal fit to similar fits using 
alternative parametrisations. All tested fit models give compatible results. The largest 
observed deviations in r sing i e are 3.9% due to modelling of the background decay time 
distribution, 0.47% due to the acceptance function and 0.39% due to the reconstructed B 
mass description, all of which are assigned as systematic uncertainties. Variations in the 
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Table 2: Corrections and systematic uncertainties on the yield ratio. 
Source Value 



Fit model 1.000 ± 0.034 

B° s mass resolution 1.000 ± 0.014 

Selection efficiency 0.968 ± 0.007 

Total correction /* 0.968 ± 0.034 



decay time resolution model are found to have negligible impact on T sing i e . 

The assumed value of the B° lifetime has a significant impact on the shape of the 
acceptance function, and the /3 parameter in particular, which in turn affects the fitted 
value of r sing i e . This effect is studied by varying the B° lifetime within its uncertainty 



25 . The largest observed deviation in r sing i e is 0.52%, which is taken as a systematic 
uncertainty. 

The fit method is tested on simulated data using large sets of pseudo-experiments, 
which have the same mass and decay time distributions as the data. Different datasets 
are generated using the fitted two-dimensional signal and background distributions, and 
Tsingie is then again fitted to these pseudo-experiments. The fit result is compared with 
the input value to search for possible biases. From the spread in the fitted values and 
the accompanying residual distributions, a small bias is found. This bias is attributed to 
the limited size of the background sample, and the resulting difficulty to constrain the 
background decay time parameters. A correction factor of 1.002 ± 0.002 is assigned to 
account for this potential bias. 

Due to the presence of a non-trivial acceptance function, the result of fitting a single 
exponential to the untagged B° s decay time distribution does mathematically not agree 
with the formal definition of the effective lifetime in Eq. ([I]), as explained in Ref. 32 



The size and sign of the difference between T sing i e and TjLgo depend on the values of t b q, 
Us, *4.Ar s , and the shape of the acceptance function. The difference is calculated with 
pseudo-experiments that sample the acceptance parameters, t b o and y s from Gaussian 
distributions related to their respective fitted and known values. Since ^.Ar s is currently 
not constrained by experiment, it is sampled uniformly from the interval [—1,1]. The aver- 
age difference between TjJ^ r0 and r single , obtained using the acceptance function affecting 

the long (downstream) K® sample, is found to be —0.001 ps (— 0.003 ps). A correction 
factor of 0.999 ±0.001 is assigned to account for this bias. 

The presence of a production asymmetry between the B® and B° s mesons could po- 
tentially alter the measured value of the effective lifetime, but even for large estimates of 
the size of this asymmetry, the effect is found to be negligible. 

Finally, the systematic uncertainties in the momentum and the decay length scale 
propagate to the effective lifetime. The size of the former contribution is evaluated by 
recomputing the decay time while varying the momenta of the final state particles within 
their uncertainty. The systematic uncertainty due to the decay length scale mainly comes 
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Table 3: Corrections and systematic uncertainties on the effective B° s — > J/ipK® lifetime. 

Source Value 



Background model 1.000 ± 0.039 

Acceptance model 1.000 ± 0.005 

Mass model 1.000 ± 0.004 

B° lifetime 1.000 ± 0.005 

Fit method 1.002 ± 0.002 
Effective lifetime definition 0.999 ± 0.001 

Total correction / c e Q ff rr 1.001 ± 0.040 



from the track-based alignment. Both effects are found to be negligible. 

The stability of the fit is verified by comparing the nominal results with those obtained 
using different fit ranges, or using only subsets of the long and downstream K® samples. 
All these tests give compatible results. 

The corrections and systematic uncertainties affecting the effective B° s — > J/ip K® life- 
time are listed in Table |3j The total systematic uncertainty is obtained by adding all the 
uncertainties in quadrature. 



7 Results and conclusion 

Using the measured ratio R = 0.0116 ± 0.0008 of B° s J/tpK^ and B° J/ipK$ 
yields, the correction factor fj? OTT = 0.968 ± 0.034, and the ratio of hadronisation frac- 
tions f s /fd = 0.256 ± 0.020 [33], the ratio of branching fractions is computed to be 

B(B ^J/ipKI) JcotT f s 

= 0.0439 ± 0.0032 (stat) ± 0.0015 (syst) ± 0.0034 {fs/fd) , 

where the quoted uncertainties are statistical, systematic, and due to the uncertainty in 
f s /f d , respectively. 



Using the known B° — > J/ip K° branching fraction 25 , the ratio of branching fractions 
can be converted into a measurement of the time-integrated B° s — > J/ipK® branching 
fraction. Taking into account the different rates of B + B~ and B°B° pair production 



at the T(4S) resonance T(B + B-)/T(B°B ) = 1.055 ± 0.025 25 , the above result is 



multiplied by the corrected value B(B° — > J/^pK°) = (8.98 ± 0.35) x 10 4 and gives 
B(B° S ^J/^K° S ) = 

[1.97 ± 0.14 (stat) ± 0.07 (syst) ± 0.15 (f s /f d ) ± 0.08 (B(B° J/ipK ))] x 10~ 5 , 

where the last uncertainty comes from the B° — > J/ipK° branching fraction. This result 
is compatible with, and more precise than, previous measurements [9,10 , and supersedes 
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the previous LHCb measurement. The branching fraction is consistent with expectations 
from [/-spin symmetry 10 . 



Using T sin gie = 1.75±0.12ps and the correction factor = 1.001±0.040, the effective 
Bg — > J/ipK® lifetime is given by 

eff reS 

T J/ipKl ~ Jcotv x r single 

= 1.75 ± 0.12 (stat)± 0.07 (syst)ps. 

This is the first measurement of this quantity. The result is compatible with the SM 
prediction given in Eq. 
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